Type I Superconductivity in YbSb 2 Single Crystals 
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We present evidence of type I superconductivity in YbSb2 single crystals, from DC and AC mag- 
netization, heat capacity and resistivity measurements. The critical temperature and critical field 
are determined to be T c w 1.3 K and H c w 55 Oe. A small Ginzburg-Landau parameter re = 0.05, 
together with typical magnetization isotherms of type I superconductors, small critical field values, 
a strong Differential Paramagnetic Effect (DPE) signal, and a field-induced change from second to 
first order phase transition, confirm the type I nature of the superconductivity in YbSb2. A possible 
second superconducting state is observed in the radiofrequency (RF) susceptibility measurements, 
with T c (2) « 0.41 K and H {2) « 430 Oe. 



PACS numbers: 74.25.-q, 74.70.Ad 
I. INTRODUCTION 



A long held empirical belief has been that type I super- 
conductors are generally elementary metals and metal- 
loids, while the majority of superconducting compounds 
exhibit type II behavior. Among the vast array of known 
binary and ternary superconductors, the number of sys- 
tems with type I superconductivity is notably limitedP^ 

YbSb2 was first claimed to be a type I superconduc- 
tor by Yamaguchi et alQ solely based on the shape of 
one M(H) isotherm at 0.4 K. Subsequently, a limited 
number of studies of the physical properties of YbSb2 
have been published. Among those, Sato et al. reported 
results of density functional theory (DFT) calculations, 
resistivity and de Haas- van Alphen measurements, which 
revealed a quasi- two-dimensional Fermi surface P Two 
other brief reports of resistivity under pressure^ and NQR 
measurements^ indicated that T c is suppressed under 
pressure p — 0.4 GPa, and that YbSb2 is likely a weakly 
coupled s-wave superconductor. Given the scarcity of 
type I superconducting compounds, and the lack of a 
thorough characterization of the magnetic and thermody- 
namic properties of YbSb2 , a detailed analysis of the su- 
perconducting ground state in this compound is needed. 
In this paper, we report results of DC and AC magnetiza- 
tion, heat capacity, resistivity and magnetic penetration 
depth, confirming the superconducting ground state with 
T c ps 1.3 K and H c s» 55 Oe. A discussion of the super- 
conducting parameters, based on the BCS and Ginzburg- 
Landau (GL) theories, is also provided. The shape of the 
M(H) curves and a second to first order phase transition 
in specific heat below T c , together with large differential 
paramagnetic effect (DPE), small critical field H c and 
GL parameter k <C 1/v2 provide strong evidence for 
the type I superconductivity in YbSb2. Moreover, the 
RF susceptibility data reveal a possible second supercon- 
ducting transition with T c (2) « 0.41 K and « 430 
Oe. 



II. EXPERIMENTAL METHODS 

YbSb2 single crystals were synthesized by flux growth 
technique, using excess amount of Sb. Elemental Yb 
(Ames Lab, 99.999%) and Sb (Alfa Aesar, 99.9999%) 
pieces in an atomic ratio of 1:9 were packed in an Alumina 
crucible and sealed in a quartz ampoule under partial Ar 
pressure. The ampoule was heated up to 650 °C, kept at 
that temperature for four hours, then slowly cooled down 
to 620 °C, after which the excess flux was removed in a 
centrifuge. The as-grown crystals were thin plates with a 
typical dimension of 5x5x0.2 mm 3 . A 1:1:1 HC1-HN0 3 - 
H 2 solution was used to remove the remnant flux from 
the surface of the crystals. 

Room temperature powder X-ray diffraction (XRD) 
measurements were carried out on a Rigaku D/Max 
diffractometer with Cu K a radiation and a graphite 
monochromator. Rietveld analysis was performed to de- 
termine the lattice parameters, using the RIETICA soft- 
ware package 

DC magnetization measurements were performed in a 
commercial Quantum Design Magnetic Properties Mea- 
surement System (QD MPMS) with a He3 insert for tem- 
peratures between 0.5 and 2 K. For the plate-like crys- 
tals, the shape was assumed to be ellipsoidal, and the de- 
magnetization factor Nd was determined from tabulated 
valuesP^l AC magnetization was measured in a dilution 
refrigerator, using a standard AC susceptometer consist- 
ing of two oppositely wound pick-up coils. An external 
modulation field of 0.1 Oe and 113.7 Hz was applied in 
the direction parallel to the crystal plate and data were 
acquired by a lock-in amplifier. After background sub- 
traction, the phase was shifted according to the excita- 
tion frequency. In order to obtain absolute values of the 
magnetization, the data were matched to the results from 
the QD MPMS measurements. The imaginary part, x" > 
was set to zero at temperatures above T c , using an ap- 
propriate offset. The offset and the scaling to absolute 
values were the same for both the temperature and field 
sweeps. 
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Heat capacity of YbSb2 was measured in a QD Phys- 
ical Property Measurement System (PPMS) with a He3 
option, using a thermal relaxation technique. To demon- 
strate the field dependence of heat capacity, measure- 
ments were carried out in a magnetic field applied per- 
pendicular to the crystal plate, with magnitudes ranging 
from to 90 kOe. Temperature dependent AC resistivity 
was also measured in the QD PPMS, utilizing a standard 
four-probe method. The sample was cut into a bar-like 
shape and four platinum wires were attached to the flat 
surface using Epo-Tek H20E silver epoxy. An AC current 
of i = 0.1 mA and / = 1000 Hz was applied along the 
in-plane direction and resistivity data were taken during 
cooling. 

In-plane magnetic penetration depth AA(T) was de- 
termined using a self-resonating tunnel-diode oscillator 
(TDR)lS operating at 16 MHz with an amplitude of 
H ac ~ 10 mOe, temperatures down to 50 mK and in 
static magnetic fields up to Hoc — 400 Oe. The sample 
was mounted on a sapphire rod with the crystal plate 
perpendicular to both Hac and Hoc- Placing the sam- 
ple into the inductor causes the shift of the resonant fre- 
quency A/(T) = — G4irx(T), where G is a calibration 
constant determined by physically pulling the sample out 
of the coil. With the characteristic sample size R, A A can 
be obtained from 4tt X = (X/R) tanh(i?/A) - lPESl 



III. RESULTS AND DISCUSSION 

The powder X-ray diffraction pattern of YbSb2 is 
shown in Fig. [T] The pattern was refined using a ZrSi2- 
type orthorhombic structure with space group Cmcm 
and lattice parameters a = 4.554 A, b = 16.715 A and 
c = 4.267 A, in good agreement with the previously re- 
ported values.^ A small amount of remnant Sb flux (~ 
5%) was found and is marked by a blue asterisk in Fig. 

m 

Fig. [2] shows the temperature dependent DC mag- 
netic susceptibility \ of YbSb2, measured in a field H 
parallel to the crystal plate for H = 10 and 20 Oe. 
The demagnetizing effect has been taken into account 
by calculating x = Xmeasured/(1 - NdXmeasured), where 
Xmeasured = M/H is the ratio of the measured magneti- 
zation M and the applied magnetic field H . The demag- 
netization factor Nd is estimated to be 0.07, if we approx- 
imate the shape of the plate-like crystal as an ellipsoid! 
For both zero field cooled (ZFC, solid lines) and field 
cooled (FC, dashed lines) data, the low temperature sus- 
ceptibility shows a clear Meissner signal at temperatures 
below T c = 1.25 K for H = 10 Oe. The superconducting 
volume fraction estimated from the ZFC data at this field 
value is very close to 100%, indicative of bulk supercon- 
ductivity. As expected, T c is suppressed by increasing 
magnetic field. The DC magnetization isotherms M(H) 
for T = 0.5 and 1.0 K are shown in Fig. [3^,, before (solid 
symbols) and after (open symbols) the demagnetization 
correction. It is clear that the corrected M(H) curves 
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FIG. 1: Powder XRD pattern of YbSb 2 (black line), with peak 
positions (vertical marks) calculated based on space group 
Cmcm and lattice parameters a — 4.554 A, b = 16.715 A and 
c = 4.267 A. A small amount of residual Sb flux is marked 
by a blue asterisk. 




FIG. 2: Temperature dependent susceptibility of YbSb2, mea- 
sured in a DC field parallel to the crystal plate. The ZFC and 
FC data are plotted as solid and dashed lines, respectively. 



show a step-like jump to zero near the critical field, char- 
acteristic of type I superconductivity. The full magneti- 
zation loops (Fig. [3b)also have the shape typical of type 
I superconductors^*^! 

The AC susceptibility x' of YbSb2 as a function of tem- 
perature is shown in Fig. [4j As the field increases from H 
— to 55 Oe, the onset temperature of the Meissner sig- 
nal drops from 1.41 K to 0.14 K. In a H = 60 Oe field, the 
superconductivity is further suppressed and cannot be 
detected down to 0.06 K. Similarly, a suppression of H c 
with T from 56.3 Oe (T = 0.06 K) to 2.3 Oe (T = 1.28 K) 
is illustrated by the x'(H) data in Fig. [5] The onset crit- 
ical temperature T c and critical field H c values from AC 
susceptibility measurements are summarized in a H — T 



phase diagram in Fig. 10 and will be discussed later. A 
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FIG. 3: (a) Isothermal magnetization and (b) full M (H) loops 
of YbSb 2 for T = 0.5 K (squares) and 1 K (triangles) in H \\ 
plate. Solid and open symbols correspond to data before (top 
axis) and after (bottom axis) the demagnetization correction, 
respectively (see text). 



noteworthy feature of the x'(T) an d x'(H) curves is the 
positive peak in the vicinity of the superconducting tran- 
sition for H > 0, known as the differential paramagnetic 
effect (DPE) in superconductors!^ The DPE signal orig- 
inates from the positive dM/dH values at temperatures 
right below H c , and can be observed in either type I or 
type II superconductors. Nevertheless, since the height of 
the DPE peak in type II superconductors cannot exceed 
the absolute value of the diamagnetic susceptibility,^ the 
observed large DPE signal (Figs. [4] and [5]) clearly points 
to type I superconductivity in YbSb2. 

The field-dependence of the superconducting transi- 
tion of YbSb2 is further confirmed by heat capacity mea- 
surements in fields up to H — 90 kOe. For clarity, a 
subset of these data are shown in Fig. [6] In the H = 
curve, a clear jump in the heat capacity confirms the 
bulk superconductivity. The transition temperature T c , 
determined at the midpoint of this jump, is close to 1.32 
K and agrees well with the previous report.^ As H in- 
creases, T c monotonically decreases and drops below 0.4 
K (the minimum available temperature for these mea- 
surements) at H = 60 Oe. The peak at the transition 
also becomes sharper and higher for H = 10 Oe, com- 
pared to that for H — 0, indicating a change from second 
to first order phase transition, commonly seen in type I 
superconductors. As the field is further increased, a non- 
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FIG. 4: Temperature dependent AC susceptibility of YbSb2, 
measured in fields H up to 60 Oe. 




FIG. 5: Field dependent AC susceptibility of YbSb 2 , 
sured for various temperatures from 0.06 K to 1.28 K. 



monotonic change of the electronic and phonon specific 
heat coefficients 7 and (3 is revealed by the C p /T vs. T 2 
plots (Fig. [6j left inset). While this field dependence re- 
mains to be understood, it makes it difficult to determine 
the electronic specific heat jump AC e / = C e i_ s — C e i^ n at 
T c . An additional difficulty in estimating AC e ; is a pos- 
sible second superconducting transition at lower temper- 
atures, which will be discussed below. We therefore use 
the H = Cp/T vs. T 2 data (black, Fig. |6j left inset) to 
estimate 7 and /? from the linear fit. This gives 7 ~ 3.18 
mJ/mol K 2 and /3 « 0.90 mJ/mol K 4 . After subtracting 
the phonon contribution C p h — /3T 3 , the jump in the spe- 
cific heat at the superconducting transition is estimated 
to be AC e i/jT c rj 0.85 (Fig. [6j right inset), smaller 
than the BCS value of 1.43. As already mentioned, this 
could be due to a second superconducting transition and 
the non-monotonic field dependence of 7 and ft, as il- 
lustrated by a subset of C P /T(T 2 ) curves shown in left 
inset, Fig. [6j The Debye temperature can also be esti- 
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FIG. 6: Temperature dependent heat capacity of YbSb2 for 
H = 0, 10, 20, 30, 40, 50, 60 Oe and 20, 60 kOe; left inset: 
Cp/T vs. T 2 for H = 0, 60 Oe, 20 kOe and 60 kOe, showing 
a field dependent Sommerfeld coefficient 7; right inset: the 
superconducting (squares) and normal (black line) electronic 
specific heat for H — 0, plotted as C e i/T vs. T. An entropy 
conservation construct (red lines) is used to determine the 
jump in the electronic specific heat at T c . 



mated, using Op, — (127r 4 7V J 4rA:B/5/3) 1/ ' 3 , where r = 3 is 
the number of atoms per formula unit. This yields Op, ss 
186 K. 

The AC resistivity of YbSb 2 , for H = and i | plate, 
is shown in Fig. [7J At high temperatures, a linear 
temperature dependence of p(T) is evident, as expected 
for metals. As seen in the upper inset, the resistivity 
drops to zero at T c — 1.25 K, with a residual resistivity 
Pq(2K) — 0.53 pVLcm just above T c . The residual re- 
sistivity ratio, calculated as RRR = p(300K) / p(2K) , is 
around 186, indicative of a good quality metal. At tem- 
peratures below 8 K and above T c , the resistivity shows 
a quadratic dependence on temperature: p = po + AT 2 . 
From a linear fit of Ap = p — po vs. T 2 , the coefficient A 
is determined to be 0.0013 pVl cmK~ 2 . The Kadowaki- 
Woods (KW) ratio A/7 2 = 12.8a , where a = 10' 5 pfl 
cm mol 2 K 2 mJ -2 is a nearly universal value observed in 
strongly correlated electron systemsP^ This large KW 
ratio is consistent with the analogous value previously 
reported 8 , and could be associated with electron-phonon 
scattering or enhanced electron correlations. 

Based on the resistivity and heat capacity data, sev- 
eral superconducting parameters, such as the London 
penetration depth A^, coherence length £, Ginzburg- 
Landau parameter n and electron-phonon coupling con- 
stant \ e i- p h can be estimated. With four formula units 
per unit cell volume (V), the electron density n of YbSb 2 
can be calculated as n = 8/V = 2.483 x 10~ 2 A~ 3 , assum- 
ing the valence of Yb to be 2+. The Fermi wave vector 
kp can be roughly estimated if we assume a spherical 
Fermi surface, which gives k F = (3wr 2 )3 = 0.903 A" 1 . 




FIG. 7: AC resistivity of YbSb 2 , for H = and i \\ plate. Top 
inset: a zoomed-in view of the low temperature resistivity 
around T c \ bottom inset: Ap = p — p vs. T 2 (symbols), 
together with a linear fit (solid line). 



The Fermi wave vector kp, together with the Sommerfeld 
coefficient 7 = 3.18 mJ/molK 2 = 6.56xl0" 5 J/cm 3 K 2 , 
yield the effective electron mass m* = h 2 k 2 F ^i j-n: 2 nk 2 B = 
1.39 mo, where mo is the free electron mass. The London 
penetration depth A^(0) and coherence length £(0) can 
also be derived as Al(0) = (to* / pone 2 ) 1 / 2 — 40 nm and 
£(0) = 0.18S 2 k F /k B T c m* = 826 nm. It results that the 
GL parameter k = A L (0)/£(0) = 0.05 < l/%/2, confirm- 
ing the type I superconductivity in YbSb2- According to 
the McMillan theory,^ the electron-phonon coupling is 
given by 

% _ 1.04+m* ln(g c /1.45T c ) 

*el-ph — (l_o.62/i*) ln(0 D /1.45T C )-1.04 

where the Coulomb pseudopotential p* is usually be- 
tween 0.1 and 0.15. Assuming p* = 0.13, the electron- 
phonon coupling is estimated to be \ e i- P h ~ 0.51, sug- 
gesting that YbSb2 is a weakly coupled BCS supercon- 
ductor. Moreover, the A e ;_ p ^ value confirms the effec- 
tive electron mass to* as calculated using m* = (1 + 
^ei-ph)iTiQ which gives to* = 1.51mo. 

Figure [8] shows frequency A/(T) ~ x(T) measured for 
a 0.7 x 0.5 x 0.3 mm 3 sample. Data above T c represent 
a combination of magnetic and resistive responses in the 
normal state of YbSb2. The skin depth at T = 1.5 K was 
estimated to be 5 ~ 8.9 pm, which is much smaller than 
any dimension of the sample. The skin depth 6 was cal- 
culated with the residual resistivity po = 0.53 pQcm and 
operating frequency of 16 MHz. A slight upturn before 
the superconducting transition can be attributed to the 
response of some paramagnetic impurities. T c was de- 
termined as the temperature where A/(T) deviates from 
the normal state behavior. 

In the pure Meissner state (Fig. [8]), for Hoc = 0, 
both ZFC and FC curves coincide. Additionally, apart 
from the superconducting transition, a small feature near 
0.11 K was observed, as shown in the inset. This may 
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FIG. 8: Shift of resonant frequency A/ = — G4n\(T), mea- 
sured in fields up to 400 Oe. The ZFC and FC data are shown 
as solid and dashed lines, respectively; inset: a zoomed-in 
view of the low temperature data for H — 0, with ZFC and 
FC data shown as solid and open symbols. 




FIG. 9: Superfluid density calculated as p = A 2 (0)/A 2 (T). 
The solid red line shows a fit to the typical behavior of an 
s-wave BCS superconductor. 



be attributed to a phase associated with extrinsic mag- 
netic impurities. In Fig. [9j the calculated superfluid 
density, p — A 2 (0)/A 2 (T), is found to be consistent with 
a single-gap s-wave pairing in YbSb2, except for the im- 
purity contribution which modifies the curve at the low- 
est temperatures. In the presence of magnetic impurities 
with magnetic permeability p(T), the actual measured 
penetration depth is renormalized as A m = \J p(T)\, so 
that the experimentally constructed superfluid density 
is given by p m — X 2 (0)/[p(T)X 2 (T)]. For paramagnetic 
impurities behaving following the Curie-law, this renor- 
malization leads to a downward trend of p m in the region 
where real superfluid density is already flat. At higher 



FIG. 10: H-T phase diagram of YbSb2. Points determined 
from x'(T), x (H), C P (T) and A/(T) data are marked by 
black squares, red triangles, blue diamonds and green circles. 
The superconductive phase boundary determined from ther- 
modynamic data is illustrated by a purple line. The possible 
new superconducting state observed in the RF magnetization 
measurement is marked with a dashed line; inset: a zoomed-in 
view of the data points from thermodynamic measurements. 



TABLE I: Measured and calculated superconducting and 
thermodynamic parameters of YbSb2. 

To (K) 



Hc(0) (Oe) 

7 (mJ/mol K 2 ) 

P (mJ/mol K 4 ) 

&d (K) 

AC el /jT c 

RRR 

A {pncm/K 2 ) 

A/>y 2 (10~ 5 ^mol 2 K 2 mJ- 

k F (A" 1 ) 

m* (mo) 

A L (0) (nm) 

€(0) (nm) 



^el — ph 

rj 2) (k) 

H { c 2) (0) (Oc) 



1.30 ± 0.2 
55 ± 2 
3.18 
0.90 
186 
0.85 
186 
0.0013 
12.8 
0.903 

1.45 ± 0.06 

40 

826 

0.05 

0.51 

0.41 

430 



temperatures the contribution of paramagnetic impuri- 
ties rapidly vanishes, since p = 1 + Anx = 1 + C/T, 
where C is the Curie constant. 

In non-zero magnetic field, ZFC-FC A/ curves (Fig. [8]) 
show hysteresis up to H — 50 Oe. Interestingly, the FC 
data indicate stronger repulsion below an intermediate 
temperature, marked with solid triangles, which system- 
atically decreases with increasing H. This crossover no 
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longer exists above 40 Oe, and ZFC data recover stronger 
repulsion. Above 70 Oe, which is greater than H c de- 
termined by various thermodynamic measurements, data 
still show diamagnetic response without hysteresis, which 
persists up to 400 Oe. It should also be noted that the AC 
susceptibility for H = 60 Oe (Fig. |4| remains slightly dia- 
magnetic after the DPE peak vanishes, consistent with 
the results displayed in Fig. [8j The origin of the dia- 
magnetism above the bulk H c is not clear, however the 
H-T phase diagram shown in Fig. [10] is reminiscent of a 
field-dependent pairing state with multiple order parame- 
ters in the heavy fermion superconductor PrOs4Sbi2 p 3 * 24 l 
More detailed measurements are required to fully under- 
stand the rich physics of this unconventional behavior, 
and to clarify whether this could be bulk or surface su- 
perconductivity. 

The relationship between T c and H c is summarized 
in the H-T phase diagram (Fig. 10). The H c values 
determined from x'(T) (squares), XyH) (triangles), C p 
(diamonds) and A/ (circles) below ~ 80 Oe can be fit 
to the expected BCS temperature dependence H C (T) = 
H C (0)[1 - (T/T c ) 2 } (solid line, Fig. [lO] inset). This gives 
H c (0) = 55 Oe and T c = 1.30 K. The possible new su- 
perconducting state inferred from the RF magnetization 
(Fig. [8| can also be described with a similar equation 

r(2) 



H?'(T) = hP(0)[1 - {T/T^') 2 } (dashed line, Fig. llOk, 



which gives #i 2) (0) = 430 Oe and = 0.41 K. 



(2) _ 



The superconducting and thermodynamic parameters 
of YbSb2 are summarized in Table[IJ Several traits of type 
I superconductors have been observed in this compound, 
including a small GL parameter k = 0.05, typical shape 
of the DC M (H) isotherms (Fig. a strong DPE sig- 
nal in the AC magnetization (Fig.Tjjl), small critical field 
values, and a change from second to first order phase 
transition induced by magnetic field and visible in spe- 
cific heat data (Fig. [6]). All these observations provide 
proof of the type I superconductivity in YbSb2. In ad- 
dition, a possible second superconducting state at lower 
temperatures is observed in RF magnetization (Fig. |8|, 
which reveals unconventional behavior, as of yet not fully 
understood. This calls for more experiments to elucidate 
the underlying physics in YbSb2- 
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